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Abstract 



The numerical solution of the recently formulated number-projected 
Hartree-Fock-Bogoliubov equations is studied in an exactly soluble cranked- 
deformed shell model Hamiltonian. It is found that the solution of these 
number-projected equations involve similar numerical effort as that of bare 
HFB. We consider that this is a significant progress in the mean- field studies 
of the quantum many-body systems. The results of the projected calculations 
are shown to be in almost complete agreement with the exact solutions of 
the model Hamiltonian. The phase transition obtained in the HFB theory as 
a function of the rotational frequency is shown to be smeared out with the 
projection. 
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The mean-field models with effective forces have been quite successful in describing the 
gross features of the quantum many-body systems. Although, the mean-field approaches 
are appropriate for systems with a very large number of particles, they have also been 
quite useful to describe the properties of finite quantum systems, for instance the atomic 
nucleus. The ground-state properties of atomic nuclei have been well described using the 
Hartree-Fock (HF) and Hartree-Fock-Bogoliubov (HFB) mean-field approaches with various 
effective interactions However, the mean- field application to a finite system suffers from 
a fundamental problem that it leads to sharp phase transitions. The phase transition is an 
artefact of the mean-field approach and is not observed in the experimental data. The phase 
transition obtained is due to the neglect of the quantal fluctuations, which become quite 
strong for finite systems. 

There are various methods in the literature to consider the quantal fluctuations on the 
mean-field solution for the finite system. One very powerful method is through the restora- 
tion of the broken-symmetries by employing the projection methods In the present work, 
we shall consider the restoration of the gauge-symmetry associated with the particle-number. 
It is known in the HFB studies that one often obtains a phase transition from the superfiuid 
paired-phase to the normal unpaired-phase. This phase transition is due to the fluctuations 
in the particle-number, the HFB wavefunction does not have a well defined particle-number. 
In most of the analysis, the particle-number fluctuations are treated in an approximate way 
by employing the Lipkin-Nogami prescription However, it has been shown that this 

approach also breaks down at high-rotational frequencies and as a matter of fact violates 
the variational principle 0. 

The exact particle number-projection can be performed by using the gradient methods 
But this approach is numerically quite involved and has been applied to only separable 
interactions with restricted model spaces PJ^]- There has been an unresolved issue whether 
HFB like equations can be obtained with the projected- energy functional. This problem 
has been recently solved |Q and it has been shown that it is possible to obtain the HFB 
equations from an arbitrary real energy functional which is completely expressible in terms of 
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the density-matrix and the pairing-tensor. The projected-energy functional can be expressed 
in terms of the density-matrix and the pairing-tensor and one obtains the HFB equations 
with modified expressions for the pair-gap and the Hartree-Fock potential The expressions 
for these quantities acquire a relatively simple form for the case of particle number-projection 



To check the applicability of the number projected-HFB (PHFB) formalism, detailed 
numerical analysis is carried out in a simple cranked-deformed shell model Hamiltonian [§. 
Although, this model cannot be used directly to study the experimental data, but it contains 
all the basic ingredients of a more realistic model. The advantage in this model is that it can 
be solved exactly and it is possible to check the accuracy of an approximate method. We 
consider that it is quite instructive to test the number-projection method in a cranking model 
as the Coriolis forces destroy the pair-correlations and the results become quite sensitive to 
the treatment of the pairing-interaction. As we shall demonstrate, the present projection 
method reproduces almost exactly the results of the shell model calculations for all the cases 
studied. 

The model Hamiltonian consists of a cranked deformed one-body term, h' and a scalar 
two-body delta- interaction |Q. The one-body term is the familiar cranked- Nilsson mean- field 
potential which takes into account of the long-range part of the nucleon-nucleon interaction. 
The residual short-range interaction is specified by the delta-interaction. The deformed shell 
model Hamiltonian employed is given by 



H = h'+V2 



2, 



(1) 



H = hdef - (^Jx + V2, 



(2) 



where, 




m 




and 



^2 = - E E^Ai^A 



(4) 
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with A^j^j^.j = (4c'J)lm and v4^^j = (A^j^^j^^Y . For the antisymmetric-normahzed two-body 



matrix-element ( Ej ), we use the delta-interaction which for a single j-shell is given by [T^ 



^ 2(2L + 1) 



1 _i 

2 2 ^ 



(5) 



where the symbol [ ] denotes the Clebsch- Gordon coefficient. We use G = g J Rf^ir^dr as 
our energy unit and the deformation energy k is related to the deformation parameter (3. 
For the case of /iii/2 shell, fi;=2.4 approximately corresponds to P = 0.23. 

It has been shown in ref. [|| that the variation of the number projected-energy functional 
results in the HFB matrix equation 



E' 



where 



H' 



^n\n2 ~^ ^nxn2 ~^ ^n-in2 ^^nin2 ^nin2 

A* c-'* F* A* _|_ \/) 

The number- projected expressions for e^^^^jF^^^j, A„^„2 and A„^„2 are given by 

e' = \jdcp y{<P) (F(0)Tr[e'p(0)] + [1 - 2ze-^^ sin #(0)]e'C(0)) 
+ h.c. 

T = \jd<P y{<P) (r(0)lTr[F(0)p(0)] + ^[1 - 2^6"^'^ sin0p(0)]F(0)C(0)) 



A 



+ h.c 
1 



+ h.c. 



2ie-'"^ sin0(:7(0)A(0)K* 



(6) 



(7) 



(8) 



(9) 



(10) 



A = -/rf</.y(</))e-2^*C(</))A(</))-(..f, 



(11) 



with 
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and 



(0)= E (0), 

n2ni 

^nsnii.'P^ r) / y ^nin2 (*^)^nin2n3n4 5 
nin2 



P(0 
y(0 



C(0)p, 

e^''^KC*{(f)) = e^"^C^{(j))K, 

e2^'^(l + p(e2*'^-l))"\ 
1 e^'^^^Metfe^'^) 



2v y'det 



Jdgx{(p)^ 



(12) 
(13) 
(14) 

(15) 
(16) 
(17) 
(18) 
(19) 

(20) 



(21) 

The quantities p and n in the above equations are the HFB density-matrix and the pairing- 
tensor, e' in (H) are the single-particle energies of the cranked-deformed potential (j^) and v 
in (p^p^) is the uncoupled antisymmetric matrix-element of the two-body delta-interaction 



The term designated by A in (10) does not appear in the ordinary HFB formalism and it 
can be immediately shown that it vanishes for the gauge-angle = 0. This term orginates 
from the variation of the pairing-energy with respect to the density- matrix. In normal HFB, 
the pairing-energy depends only on the pairing-tensor, but the PHFB pairing-energy also 
depends on the density-matrix through the norm-overlap. Actually, in the general the norm- 
overlap depends on both density- matrix and the pairing-tensor [§. But for the special case 
of number-projection, the term in the overlap-matrix which depends on the pairing-tensor 
can be rewritten in terms of the density-matrix by using the HFB relation (p — = kh)). 
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Due to this transformation, the expression for A in ([TT|) has a very simple appearance and 
reduces to the famihar form in the canonical representation P,p]l|. 

The integration in (p|-p!l|) over the gauge-angle has been performed using the Gauss- 



Chebyshev quadrature method |T^. In this method, the integration over the gauge-angle is 



replaced by a summation. It can be shown |T2| that the optimal number of mesh-points in 
the summation which eliminates all the components having undesired particle numbers is 
given by 

M = max (-N, Q~-) +1, (22) 



.2 2. 

where N is the number of particles and Q is the degeneracy of the single-j shell. In the 
present study with N=6 and Q = 6, the optimal number of points required is M = 4. 

In the present analysis of a single-j shell, the basis in which the HFB matrix is constructed 

are the magnetic sub-states of j = 11/2 with m = (11/2,9/2, ,—9/2,-11/2). The 

summation indices ni,n2, and in the all the expressions given above run over these 
magnetic states. In order to check the dependence of the HFB and PHFB results on the 
pairing interaction, the calculations have been performed with monopole (L=0), monopole 
plus quadrupole (L=0 and 2), and with full delta- interaction. The results of the HFB and 
PHFB will be compared with the exact results for the three pairing-interactions. 

The results of the cranking calculations with monopole-pairing force are compared in Fig. 
1. In the three-pannels of the figure, we compare the total-energy (Efot), the pairing-energy 
(Epair) and the alignment (< >) which is the expectation value of the angular-momentum 
along the rotational function of the rotational frequency. The expressions for 

the total energy is given by 

Etot = J d<j)y{<j)) {H,p{<f)) + + Hp^m , (23) 

where 

i/,p(0) = Tr(ep(0)), (24) 

i7,,(0) = ^Tr(r(0)p(0)), (25) 



^Pp(0) = -2Tr(A(0)7^*(0)). (26) 
The expressions for the pairing-energy and the ahgnment are given by 

Epair = J d(l)Hpp{(l)), (27) 

<J,> ^ j d(tn^{J,p{4>)). (28) 

It can be easily shown that for the gauge-angle, = 0, the normal HFB expressions for 
these quantities are recovered. 

It is apparent from the top pannel of Fig. 1 that the results of the exact shell model and 
the PHFB are very similar, the two curves are almost indistinguishable for all the frequency 
points. The results of HFB on the other hand deviate considerably from the exact results 
at lower frequencies. However, it can be seen from Fig. 2 that the HFB results converge 
towards the exact results with increasing rotational frequency. The HFB energy before the 
bandcrossing at Tiuj — 0.5G is shifted from the exact energy by a constant amount and can 
be improved by renormalising the strength of the pairing-interaction. Therefore for the total 
energy, the HFB approach is not a poor approximation. The actual problem in HFB lies 
in the analysis of the pairing-energy and the alignment which are shown in the two lower- 
pannels of Fig. 1. The HFB pairing-energy has a finite value till Tuo — 0.45G and then 
suddenly goes to zero at fiu) = 0.5G. This transition is an artefact of the HFB approach 
as is clearly evident from Fig. 1. The PHFB pairing-energy does drop at the bandcrossing 
but has a finite value at all the rotational frequencies. This feature of the pairing- energy is 
also reflected in the alignment, which is determined by the competition between the Coriolis 
forces and the pairing- correlations. The alignment till fiuj — 0.45G is very similar in all the 
three cases but then HFB value suddenly jumps to < >— 10, the exact and the PHFB 
values on the on other hand do not show this sharp transition. The gain in alignment at 
the first-crossing, referred to as the (AB)-crossing, is quite similar in all the cases, however, 
in HFB the crossing occurs somewhat earlier. 

The results of pairing-interaction with monopole plus quadrupolc terms are shown in 
Fig. 2. The comparison among HFB, PHFB and the exact results is quite similar to Fig. 2. 
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The total HFB energy is shifted by a constant factor from exact and PHFB energies before 
the bandcrossing. After the crossing, the HFB results become closer to the exact. The 
HFB pairing-energy depicts a transition and PHFB pairing-energy on the other hand drops 
smoothly at the bandcrossing. However, the HFB gain in ahgnment at the (AB)-crossing is 
much lower than 10 and has clearly a wrong behaviour after the crossing. 

The results with the full delta-interaction are presented in Fig. 3. The HFB total-energy 
is now in a better agreement with exact and PHFB as compared to the results shown in Figs. 
1 and 2. In fact, it is evident by comparing the three figures that the total-energy improves 
by including higher multipoles in the pairing-interaction, maximum deviation is noted for 
the monopole case. This appears to be in contradiction to our basic understanding of the 
mean-field approach, one normally expects that HFB or BCS is a better approximation for 
the pure monopole-pairing. However, it should be mentioned that in our HFB and PHFB 
analysis, we have included all the terms in the Hamiltonian. In particular, the particle-hole 
contribution (P) amounts to about 6G in the total-energy and is maximum with the full 
delta-interaction. If one excludes this contribution, the discrepancy would be largest for the 
delta- interaction. 

The pairing-energy in Fig. 3 again depicts a phase transition at hcu — 0.55G which is 
slightly higher than with monopole-interaction. The HFB (AB)-crossing with full delta- 
interaction is now close to the exact and PHFB. The gain in alignment at the (AB)-crossing 
is lower than 10 as in the case of monopole and quadrupole pairing force. The overall 
agreement with full delta-interaction appears to be better for HFB. The good agreement 
between PHFB and exact on the other hand is independent of the interaction used. 

To conclude, in the present work the number-projected HFB approach recently developed 
has been applied to an exactly soluble cranked-deformed shell Hamiltonian. The main 
motivation has been to check the numerical applicability of the projection method. It is 
clear from the present study that the projected Hartree-Fock-Bogohubov approach gives an 
accurate description of the yrast-statcs of the model Hamiltonian. The transition from a 
superfluid to the normal phase obtained in the HFB theory is shown to be smeared out with 
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the projection. 

We would like to stress that the main advantage of the present projection method is 
that it has the same structure as that of normal HFB equations. Therefore, one can use the 
existing HFB computer codes and only the expressions for the Hartree-Fock potential and the 
pairing-field need to be redefined. Instead of the normal HFB fields, in the projection method 
one needs to calculate the projected quantities as given by (|8|-|TI1). In the present model study, 
we find that the numerical work involved in the projection is similar to performing the 
bare HFB calculations. For each rotational frequency, the average CPU time on Pentium 
(166MHz) was 6.14s with projection as compared to 5.97s for normal HFB. The present 
projection method, therefore, preserves all the mathematical and computational simplicity 
of the HFB mean-field approach. 



9 



REFERENCES 

[1] Ring P. and Schuck P., The Nuclear Many-Body Problem (Springer Verlag, New York, 
1980) 

[2] Magierski P., Cwiok S., Dobaczewski J., and Nazarewicz W., Phys. Rev. C 1993. V. 48. 
P. 1681. 

[3] Gall G., Bonche B., Dobaczewski J., Flocard H., Heenen P.-H., Z. Phys. A 1994 V. 348. 
P. 183. 

[4] Valor A., Egido J.L., and Robledo R., Phys. Rev. C 1996. V. 53. P. 172 
[5] Flocard H. and Onishi N., Ann. Phys. (N.Y.) 1996. V. 254 P. 275 
[6] Egido J.L. and Ring P., Nucl. Phys. A 1982. V. 383. P. 189 
[7] Egido J.L. and Ring P., Nucl. Phys. A 1982. V. 388. P. 19 
[8] Sheikh J.A. and Ring P., Nucl. Phys. A 2000 V. 665. P. 71. 

[9] Sheikh J.A., Nagarajan N.A., Rowley N. and Pal K.F., Phys. Lett. B 1989. V. 223. P. 
1. 

[10] Brussaard J.P and Glaudemans P.W.M, Shell Model Applications in Nuclear Spec- 
troscopy (North- Holland, Amsterdam, 1977) 

[11] Dietrich K., Mang H.J. and Pradal J.H., Phys. Rev. 1964. V. 135. P. B22 

[12] Hara K., Iwasaki S. and Tanabe K., Nucl. Phys. A 1979. V. 332. P. 69. 



10 



FIGURES 

FIG. 1. The results of the total energy (Etot), the pair-energy {Epair) and the alignment (Jj) 
for six-particles in a deformed j = 11/2 orbitial using the monopole-interaction. The PHFB results 
are indistinguisable from the exact shell models results. 

FIG. 2. The results of the total energy {Etot), the pair-energy {Epair) and the alignment {Jx) 
for six-particles in a deformed j = 11/2 orbitial using the monopole plus quadrupole interaction. 

FIG. 3. The results of the total energy {Etot), the pair-energy {Epair) and the alignment {Jx) 
for six-particles in a deformed j = 11/2 orbitial using the full delta- interaction. 
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